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Abstract 

The exclusive production of pairs of vector mesons with J PC = 1 in e + e~ collisions can 
proceed through e + e~ annihilation into two virtual photons. At energies much greater than the 
meson masses, the cross section is dominated by the independent fragmentation of the virtual 
photons into the vector mesons. The fragmentation approximation is used to calculate the cross 
sections and angular distributions for pairs of vector mesons that can be produced at the B factories. 
The predicted cross sections for p°+p° and p°+4> production agree with recent measurements by the 
BaBar Collaboration. For the production of two charmonium vector mesons, the nonfragmentation 
corrections to the cross sections are calculated by using the NRQCD factorization formalism. The 
predicted cross sections for J/ip + J/ijj and J/tp + ip(2S) production are compatible with upper 
limits set by the Belle Collaboration. 

PACS numbers: 13.60.Le,12.40.Vv,13.66.Bc,14.40.Gx 



I. INTRODUCTION 



The production of hadrons in e + e~ collisions at energies well below the mass of the Z° 
proceeds predominantly through the annihilation of e + e~ into a virtual photon. This process 
can only produce final states with charge conjugation quantum number C = — 1. Final states 
with C = +1 can be produced by the annihilation of e + e~ into two virtual photons, but 
the cross sections for these processes are suppressed by a factor of a 2 , where a ~ 1/137 
is the QED fine structure constant. Final states with C = +1 can also be produced by 
the annihilation of e + e~ into a virtual Z°, but the cross sections for these processes are 
suppressed by a factor of (s/M§) 2 if the e + e~ center-of-mass energy ^Js is small compared 
to Mz- Thus, unless there are enhancement factors to compensate for the factors of a 2 or 
(s/Mf) 2 , the production cross sections at I?-factory energies for C — +1 states are orders 
of magnitude smaller than those for C = — 1 states. 

In e + e~ collisions with large center-of-mass energy s/s, factors of y/s/m, where m is 
a hadronic mass scale, can compensate for the suppression factor of a 2 . An example of 
such compensation can be seen in the exclusive double-charmonium production process 
e + e~ — > J/ip + J/ip U E|. At the beam energy of the B factories, \fsj1 = 5.29 GeV, the 
leading-order prediction for the production cross section for the C = +1 final state J/ip+J/ip 
has roughly the same order of magnitude as that for the corresponding C = — 1 final state 
J/ip + f] c |3j. The production cross section for J/ip + i] c scales as a 2 a 2 rrft c / 's 4 , where m c is the 
charm-quark mass. In contrast, the production cross section for J/ip + J/ip scales as a 4 /s 
because there is a photon-fragmentation contribution that corresponds to the annihilation 
process e + e~ — > 7* + 7*, followed by independent fragmentation processes 7* — > J/ip. The 
enhancement factor of (s/m 2 ) 3 in the production cross section for J/ip + J/ip relative to that 
for J/ip + r] c partially compensates for the suppression factor of (a/a s ) 2 . 

A similar enhancement mechanism is present in the exclusive production of a pair of light 
vector mesons V1 + V2, both of which have quantum numbers J PC = 1 . The production 
process proceeds predominantly through the annihilation process e + e~ — > 7* + 7*, followed 
by the independent fragmentation processes 7* — ► V{. We can compare the production cross 
section for V\ + V2 to that for V + P, where P is a light pseudoscalar meson. In the ratio 
of the former cross section to the latter, there is an enhancement factor of (s/Aq CD ) 3 that 
compensates for the suppression factor of (a/a s ) 2 . In this paper, we calculate the cross 
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sections and angular distributions for the exclusive production of pairs of vector mesons 
with J PC = 1— . 

The remainder of this paper is organized as follows. In Sec. |H1 we discuss the various 
contributions to the amplitudes for the production of two vector mesons and estimate their 
relative sizes. In Sec. Illl[ we give expressions for the differential cross sections for the pro- 
duction of two vector mesons in the fragmentation approximation and give numerical results 
for the cross sections. Sec. IIVI contains a computation of the nonfragmentation corrections 
to the cross sections for the production of two charmonium vector mesons. In Sec. 13 we 
compare our results with experimental measurements and with previous theoretical calcu- 
lations. The Appendix contains the expressions for the nonfragmentation corrections to the 
differential cross sections for two charmonium vector mesons. 

II. PRODUCTION AMPLITUDES 

The lowest-order QED diagrams that correspond to the creation of the quarks in the two 
vector mesons V\ and V2 are shown in Figs. ^ and El If the two vector mesons have different 
quark contents qq and q'q', then only the diagrams in Fig.^are present. Once the quarks are 
created, the formation of the vector mesons involves the exchange of arbitrarily many gluons 
between the quark lines. Gluons can be exchanged between the quark and antiquark that 
form a given meson or between quarks or antiquarks that form different mesons. However, 
factorization theorems and the scaling of various contributions yield a simplification of this 
picture in the limit i?beam — * 00. In what follows, we take the beam energy E heam = y/s/2 
to be much greater than the masses m v . of the vector mesons, which implies that -Ebeam is 
also much greater than Aqcd- 

Let us first consider the exchange of gluons only between the q and q within a given 
meson. Then, the diagrams in Fig. El are suppressed in comparison to those in Fig. ^ In 
the diagrams of Fig. Q the virtual-photon propagators are 1/rriy., and the two hadronic 
electromagnetic currents each give a factor of order Aqcd for light mesons and m c for char- 
monium mesons. In the diagrams of Fig. El the virtual-photon propagators are of order 
l/i^eam' an d the two hadronic electromagnetic currents each give factors of order -Ebeam- 
This leads to the following suppression factors for the diagrams of Fig. El relative to those of 
Fig. HI AQ CD /i?^ eam for two light mesons and ml/E^ eam for two charmonium mesons. We 
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have taken the light-meson masses to be of order Aqcd and the charmonium masses to be 
of order m c . (The diagrams of Fig. El do not contribute in the case of a light meson and a 
charmonium meson.) 

Next, let us consider the exchange of gluons between quarks or antiquarks in different 
mesons, with some of those gluons soft or, in the case of a light meson, collinear. In the case 
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to 



of two light mesons, one can use standard methods for proving factorization theorems 
show that the soft contributions cancel up to corrections of relative order AQ CD /Z<^ eam and 
the collinear contributions cancel up to corrections of relative order Aq CT) / E^^. 1 Factor- 
ization theorems are not well established for the production of heavy- quarkonium mesons. 
However, on the basis of existing factorization technology, it seems plausible, in the case 
of one light meson and one charmonium meson, that soft contributions cancel up to terms 
of relative order AQ CD (m c t> ) 2 /E^ eam and that collinear contributions cancel up to terms of 
relative order Aq CD /i?jf; eam . Here v is the typical velocity of the charm quark or antiquark 
in the charmonium rest frame, (v 2 ~ 0.3.) In the case of two charmonium mesons, there 
are no collinear singularities. It is plausible that soft contributions cancel up to terms of 
relative order (m c v) A / E^ eaxa in this case. 

Finally, let us consider the exchange of hard gluons between quarks or antiquarks in 
different mesons. (Hard gluons have energies and momenta of order .Ebeam-) In the case 
of the diagrams of Fig. |21 exchanges of hard gluons are suppressed as a s (E heam ). In the 
case of the diagrams of Fig. ^ it is necessary to exchange at least two gluons between 
the mesons in order to keep the qq pair in a color-singlet state. Hence, such exchanges 
are suppressed as a 2 (E\> eSum ). They are also suppressed by a kinematic factor that arises 
as follows. In the diagrams of Fig. without hard-gluon exchange, the virtual-photon 
propagators are 1/rriy. and the two hadronic electromagnetic currents each give a factor of 
order Aqcd for light mesons and m c for charmonium mesons. In the diagrams of Fig. ^ 
with hard-gluon exchanges, the virtual-photon propagators are of order 1/E^ eam , and the 



1 The results that we quote here apply to the case in which one sums over all polarizations of the final- 
state mesons. If one observes the polarization of one meson, then the soft contributions cancel up to 
terms of relative order Aq CD /ii^ eam , and the collinear contributions cancel up to terms of relative order 
Aqcd / ^bcam- If one observes the polarization of both mesons, then the soft contributions cancel up to 
corrections of relative order Aq CD /^ eam and the collinear contributions cancel up to corrections of relative 
order Aqcd / Bbcam- Similar reductions of the suppression factors for soft and collinear contributions occur 
in the cases in which one or more charmonium mesons are produced. 
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FIG. 2: QED nonfragmentation diagrams that contribute to exclusive two- vector- meson production 
in e + e~ collisions if the two vector mesons have the same quark content qq. 

two hadronic electromagnetic currents each give factors of order -Ebeam- This leads to the 
following additional suppression factors for hard-gluon exchanges in the diagrams of Fig. d 
AQ CD /Z^ eam for two light mesons, m c AQ C D/-Ebeam f° r a light meson and a charmonium 
meson, and m 2 c / 'E^ eam for two charmonium mesons. Again, we have taken the light-meson 
masses to be of order Aqcd and the charmonium masses to be of order m c . 

Let us now summarize the results of the preceding analyses. We conclude that we need 
only consider the fragmentation diagrams in Fig.^ up to terms of relative order Aq CD /i?^ eam 
for two light mesons and ml/E^ eam for two charmonium mesons. (The diagrams of Fig. 121 
do not contribute in the case of a light meson and a charmonium meson.) We also conclude 
that we need only consider the QCD corrections to the diagrams of Fig. ^ that correspond 
to the exchange of gluons between the q and q in each meson, up to terms of relative 
order A| CD /£b eain for two !ig ht mesons, A 2 QCB /El cam or a 2 s (E heam ) (m c A QCD /£2 cam ) for a 
light meson and a charmonium meson, and (m c t>) 4 /i?^ eam or a^(-E'beam) ( m c/-^bcam) f° r ^ wo 



5 



TABLE I: Masses, electronic decay widths, and conpling constants SVl for vector mesons V. All 
of the data except for those of p° are taken from Ref. |9j. The data for p° are taken from Ref. [10| 



in order to maintain consistency between the narrow-width approximation and the nonzero-width 
parametrization. The uncertainties shown for gy 7 are those that arise from the uncertainties in 
the electronic widths of the vector mesons. 



V 


Mass (MeV) 


Width (keV) 


g Vl (GeV 2 ) 


p° 


775.65 ± 0.64 ± 0.50 


7.06 ±0.11 ±0.05 


0.122 ±0.001 


CO 


782.65 ±0.12 


0.60 ±0.02 


0.036 ± 0.001 




1019.46 ± 0.019 


1.27 ±0.04 


0.078 ± 0.001 


J/V 


3096.916 ±0.011 


5.55 ±0.14 


0.860 ±0.011 




3686.093 ± 0.034 


2.48 ±0.06 


0.746 ± 0.009 



charmonium mesons. 2 

The QCD corrections that involve the exchange of gluons between the q and q in each 
meson are precisely those that correspond to the formation of a vector meson V from a qq 
pair that has been created at a point. They can be expressed in terms of a QCD matrix 
element of the electromagnetic current: 

J M (x) = e u u(x)^u(x) + e d d(x)^d(x) + e s s(x)^s(x) + e c c(x)7 M c(x), (1) 

where e u = e c = +| and = e s = — ~. The matrix element can be used to define a 
vector- meson-photon coupling constant gy 7 : 

(V(\)\r(x = 0)\0)=g v ,e»(\y, (2) 

where A is the helicity of the vector meson. The coupling constant gy 7 can be determined 

2 There are several recent papers in which the rate for the process 7*7* — * p°p° is calculated in the high- 
energy and diffractive limit s ^> — t, where s and t are the Mandelstam variables 0, S B 0- I n thl s 
limit, the dominant contributions to this process involve the exchange of gluons between the quarks and 
antiquarks that form different mesons. There are analogous contributions to the process e + e~ — > p°p°, 
since it also proceeds through 7*7* — > p°p°. However, as we have discussed above, these contributions 
are suppressed, and the dominant mechanism in e + e~ — > p°p° is independent photon fragmentation. 
Independent photon fragmentation is not considered in Refs. g g fl ff because it corresponds in the 
process 7*7* — > p°p° to a disconnected diagram. 
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from the electronic width of the vector meson: 

T[V - e V] = — -fl. (3) 

Using the measured electronic widths that are given in Ref . \9] , we obtain the values of gy 1 
that are shown in Table U] The stated uncertainties arise from the electronic widths. The 
uncertainties in the masses are at most 0.06% and can be neglected. 



III. CROSS SECTIONS IN THE FRAGMENTATION APPROXIMATION 

We express our results for the cross section for e + e~ — ► V\ + V 2 in terms of dimensionless 
variables ry 1 and ry 2 , which are defined by 

If we take y/s to be 10.58 GeV, the center-of-mass energy of the B factories, then r\ is small 
for charmonium {t 2 j,^ = 0.086) and very small for light vector mesons (r 2 = 0.0054). The 
differential cross sections da jd cos 9 for each of the helicity states of V\ and V 2 are 

£[H(AJ + V 2 (A 2 )] = ^^^'"M^MruY) ( 
<fa s 5 r;^[(l-x 2 )A(l,rf,rl)+4r;^] 

where x = COS6 1 , = r^, and 

A(a, 6, c) = a 2 + 6 2 + c 2 - 2ab - 26c - 2ca. (6) 

The functions i 7 \ 1> A 2 ( r i) r 2, x), which depend on the helicities of the vector mesons, are given 
by 

^±i,Ti( r i> r 2, x) = f+-(r 1 , r 2 , x) + /+-(ri, r 2 , -x), (7a) 

F±i )0 (r 1 , r 2 , x) = f + o(n, r 2 , x) + / + o(n, r 2 , -x), (7b) 

^o,±i(n,r 2 ,x) = / +0 (r 2 ,ri,x) + / +0 (r 2 ,ri, -x), (7c) 

^±i,±i(n,r 2 ,x) = / ++ (ri,r 2 ,x) + f ++ (n,r 2 , -x), (7d) 

^o,o (ri,r 2 ,x) = 16r 2 r 2 x 2 (l -x 2 ), (7e) 



where 



/ + _(n,r 2 ,x) = |(l + x)(l-x) 3 (l-r 2 -r 2 ) 2 , (8a) 
/+o(ri, r 2) x) = r 2 (l - x) 2 [(l - r 2 )(l + x) - r 2 (l - x)] 2 , (8b) 
f ++ (n, r 2 , x) = i(l - x 2 )[(l + x)r 2 (l - r 2 ) - (1 - x)r 2 (l - r 2 ) + 2r 2 r 2 x] 2 . (8c) 



Since t\ and r 2 are small, the cross section is largest if the two vector mesons are transversely 
polarized with opposite helicities: Ai = — A2 = ±1. If the vector meson Vi is longitudinally 
polarized, then the cross section is suppressed as rf. If the vector mesons have the same 
helicities, then the cross section is suppressed as four powers of r\ or r 2 . 

After summing over the helicities Ai and A2 of the vector mesons, we obtain 



da, , 32 7 r 3 aX 7 ^ 27 A 1/2 (l, r?,r* 

— {m Vl ,m V i - 



16r 2 r 2 (r 2 + r\ 



dx K Vli s 5 [ ri r 2 (l -r\ -rl)] 4 [l - (1 - A)x 2 ] 2 

+2(1 - * 2 )(1 + rf + r 2 4 )A(l, r 2 , r 2 ) - (1 - * 2 ) 2 A 2 (1, r 2 , r 2 ) 



(9) 



where A = 4r 2 r|/ (1 — r 2 — r|) 2 . The numerical value of A is very small for both charmonia 
and light vector mesons. For example, A = 0.043 for J/ijj + J/ijj and A = 0.00012 for 
p° + p°. In the arguments on the left side of Eq. 0, we have indicated the dependences on 
the masses of the vector mesons, as they will play a role in a subsequent discussion of the 
effect of the width of the p meson. The total cross section a is obtained by integrating over x 
from —1 to +1 if the two vector mesons are distinct. If the two vector mesons are identical, 
then the integration range of x is taken to be to 1 in order to avoid double counting. The 
integration formulas that are required to evaluate the definite integrals over x from —a to a 
analytically are 



, [1 - (1 - A)x 2 ] 2 1 - a 2 + a 2 A + 2v ^A 1 - av^A ' 1 J 
1 - x 2 a A 

dx-, 



[1 - (1 - A)a; 2 ] 2 (1 - A)(l - a 2 + a 2 A) 



2- A , 1 + aVl- A 
+ 2(l-A)3/ 2 ln i- aV l-A ' (1 ° b) 



(1-x 2 ) 2 1 / aA 2 
dx—^- +— - = — 2a + 



[1 - (1 - A).x 2 ] 2 (1-A) 2 V l-a 2 + a 2 A 

+ A(A :ji* 1+av ^- (ioc) 

2(l-A) 5 /2 i_ av ^A V ; 

The cross sections for two vector mesons with J PC = 1 have also been calculated 
recently in Ref. jll| . As was pointed out in Ref. |11], at the level of precision of this work, 
it is necessary to take into account the nonzero width of the p meson. Following Ref. jll|. 
we note that Eq. Q can be generalized to the case of a continuous mass spectrum: 



r/cj = \4n 2 a) J dm i dm 2 a e+e--+H(mi) a e + e -^ H (rn 2 ) ^ 2 (mi, m 2 ), (11) 



where <r e + e -_ > #(m) is the total cross section for the process e + e~ — > hadrons at an energy 
m in the e + e~ rest frame, and da 7 * 7 */dx is the cross section for e + e~ annihilation into two 
massive photons, da y * y */dx can be related to da/dx in Eq. (JOJ): 



da. 



/ ~ 2 ~ 2 



dx 



' 2 'mi,m 2 ) = ( - — "" L "" z ) rL( mii m 2 y 

\ / ±'Kagv 11 gv2i J dx 



da 



(12) 



The ranges of integration of rh\ and m 2 in Eq. (|TT|) are determined by the physical region 
for the two-meson final state: fhi < y/s — my i for the case of a p° meson and another meson 
Vi and rh\ + m 2 < for the case of two p° mesons. One can recover the narrow-width 
approximation by writing 



a 



NW 



/ 4irag Vr 



7r5(m —m v ). 



(13) 



We parametrize the contribution of the p° meson to a e + e -^ H (m), following Ref. [121 ]. as 

8vra 2 



o- e +e-^ P {m) 



-^pl(m 2 )e(m - 2m 7r )|F 7r (m 2 )| 2 , 



where is the pion mass and p n (s) is the pion momentum: 

p^s) = a/s/4 - m 2 . 
The form factor ^(s), which depends on the parameters M p , T p , and /3, is 



(14) 



(15) 



FJs) 



Ml-s + f(s) - iM p T p {s) 



where 



d 



3ml M p + 2^(M 2 ) 
— log p — + 



Trpl(M? 



2m„ 



2vr^(M 2 ) 7rpl(M ; 



mlM p 



(16) 



(17) 



r„(s) is the energy-dependent width of the p, 



r P 00 = r, 



LrtMJ 



M 2 



1/2 



(18) 



and the real part of the correction to the denominator is given by 

Ml 



m 

h(s) 
h'(s) 



{pl(s)[h(s) - MM 2 )] + (M 2 - s)pl{Ml)h\Ml)} 



— log 

2m,r 



7T V s 

dh(s) 
ds 



(19a) 
(19b) 
(19c) 



TABLE II: Cross sections in units of fb for e + e~ — ► V\ + Vi at -Ebeam = 5.29 GeV, calculated 
by using the fragmentation approximation. The uncertainties that are shown are only those that 
arise from the uncertainties in the electronic widths of the vector mesons. The first five rows are 
calculated in the narrow-width approximation. The last row is calculated by taking into account 
the nonzero width of the p meson, as is described in the text. 



Vi \ v 2 


P° 




4> 




if>(2S) 


p° 


139.61 ±4.82 


23.47 ± 0.88 


35.93 ± 1.29 


41.67 ± 1.27 


15.60 ±0.46 






0.99 ± 0.07 


3.02 ±0.14 


3.50 ±0.15 


1.31 ±0.05 








2.30 ±0.15 


5.20 ±0.21 


1.94 ±0.08 


J/v 








2.52 ±0.13 


1.81 ±0.06 












0.32 ±0.02 


P° 


126.08 ±4.36 


22.31 ±0.84 


34.18 ± 1.23 


39.83 ± 1.22 


14.92 ± 0.44 



The parameters M p and V p are obtained from a fit to the data for the e + e — > 7t + tt cross 
section Q: M p = 775.65 ± 0.64 ± 0.50 MeV and T p = 143.85 ± 1.33 ± 0.80 MeV. (Similar 



results have been obtained in Ref. 



13j . ) The parameter (5 is not given in Ref. 
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. We infer 



it from the value T[p — > e + e ] = 7.06 ± 0.11 ± 0.05 MeV that is given in Ref. 10] and the 
formula for the electronic width that is given in Ref. |l2|: 

^ e+e -i _ ^ 2 pl(K)d + dr P /M p y 
1 [p e e 1 - 9m p t p (i + ■ (2Uj 

Using the values for M p and T p from Ref. |l0|, we obtain j3 = —0.0815234. 

In Table ILTj we give our results for the integrated cross sections, calculated in the narrow- 
width approximation, for the production of V1 + V2 for the vector mesons Vi = p°, u, <fi, J/ip, 
and ip(2S). We also give integrated cross sections for the production of p° + V in which the 
nonzero width of the p has been taken into account. The uncertainties that are shown are 
only those that arise from the uncertainties in the V — 7 coupling constants that are given in 
Table |1] The effect of the nonzero width of the p meson is to decrease the cross sections for 
p° + p° by about 10%, for p° + u and p° + by about 5%, and for p° + J/ip and p° + ip(2S) 
by about 4%. 

We note that the differences between the zero-width and nonzero-width results arise 
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FIG. 3: Differential cross sections do j d\x\ for e~^~e annihilation into V\ + V2 at £<beam — 5.29 GeV, 
where VI = p , <j), or J ftp. The areas under the curves are the total cross sections. 



mostly from the fact that the quantity 



dm 2 a, 



m 



irrii 



(21) 



differs from the coefficient of the 8 function in the narrow-width approximation of Eq. (|13j). 
The quantity I p is about 6% smaller than the coefficient of the 8 function in Eq. (|13j). (The 
effect of the nonzero width of the p meson itself is to actually increase the cross sections.) 
The coefficient of the 8 function in Eq. (1 1311 derives from the electronic width of the p 
meson, which, in turn, is calculated in Ref. [lj| by using Eq. (120(1 . Eq. (f2"UJ) is derived from 
the vector-meson-dominance model Q|. Hence, the experimental electronic width of the p 
meson that is given in Ref. depends on the assumptions of that model. An alternative 
definition of the electronic width of the p meson can be obtained by equating I p to the 
coefficient of the 8 function in Eq. (fTTjJ) and using Eq. (j3J) to relate gy a to Y[p — > e + e~]. This 
approach leads to the result 



m. 



e e 



oh- (22) 

differs from the one in Eq. (I20j) by terms of 



It can be shown that the definition in Eq. 
order Y 2 p /m 2 p . 

In Fig. 01 we show the differential cross sections da/d\x\ for the production of the identical 
vector mesons p° + p°, 4> + (ft, and J/ip + J/ip and for the production of the distinct vector 
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mesons p° + 0, p° + J/ip, and <p + J/ip. The differential cross sections peak sharply near 
| a; | = 1 before falling to zero at the endpoint. For example, the maximum values of the 
differential cross sections are 7.1 x 10 4 fb at x = 0.99994 for p° + p° and 13.3 fb at x = 0.993 
for J/V + J/ip. 

We now discuss the theoretical uncertainties in the predictions for the cross sections in 
Table |H] The scaling of the theoretical uncertainties was determined in Section |Hj In the 
case of two light vector mesons, the leading theoretical uncertainties arise from the fragmen- 
tation approximation and from exchanges of collinear gluons between the mesons. These 
uncertainties both scale as (AQCD/^beam) 2 - If we take Aqcd ~ 0.5 GeV, then the estimated 
fractional uncertainty is 0.9%. In the case of a light vector meson and a charmonium meson, 
there is an additional theoretical uncertainty that arises from the exchange of hard gluons 
between the mesons that scales as a 2 (E^ eam ) (m c AQCD/-Ebeam)- If we take a 2 (E^ eam ) = 0.25 
and m c = 1.4 GeV, then the estimated fractional uncertainty is 0.2%. In the case of two 
charmonium mesons, the leading theoretical uncertainty arises from the fragmentation ap- 
proximation and scales as m 2 c / E^ eSirn . If we take m c = 1.4 GeV, then the estimated fractional 
uncertainty is 7%. This uncertainty can be reduced by calculating the contributions to the 
cross section from nonfragmentation diagrams. We discuss this calculation in Sec. IIV1 

IV. NONFRAGMENTATION CORRECTIONS TO THE PRODUCTION OF 
TWO CHARMONIA 

In this section, we calculate corrections to the production cross sections for two charmo- 
nium mesons that arise from the nonfragmentation contributions to the amplitudes. The 
contributions from the nonfragmentation diagrams in Fig. El can be calculated using the 
NRQCD factorization formalism This formalism allows QCD radiative corrections and 
relativistic corrections to be taken into account systematically. For example, the expression 
for the V — 7 coupling constant defined by Eq. (jSJ), including the leading QCD radiative 
correction and the leading relativistic correction, is 



Here if)i and x are the two-component Pauli operators in the NRQCD formalism that create 
a heavy quark and antiquark, respectively. (v 2 )v is proportional to a ratio of matrix elements 




(23) 
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of NRQCD operators between the state V and the vacuum. The factor e c comes from the 
electromagnetic current in Eq. ([T}. The factor \J1my takes into account the difference 
between the standard relativistic and nonrelativistic normalizations of the state |V(A)). 

The NRQCD factorization formalism was used in Refs. to calculate the cross sections 
for double-charmonium production from e + e~ annihilation into two virtual photons. The 
calculation of the cross sections and the estimates of theoretical errors were carried out in a 



way that was as close as possible to a previous calculation of the cross sections for double- 

f] 

charmonium production from e + e annihilation into a single virtual photon 3]. The cross 
sections were expressed in terms of the coupling constants a and a s , the pole mass m c of the 
charm quark, and a factor (0\)y that is related to the NRQCD matrix element in Eq. (|23jl : 

(Oi)v = lY,\( V ( X ^x\°)\ 2 - ( 24 ) 

A 

The fragmentation terms in the cross sections in Ref. Q] can be recovered by replacing o£ 
in Eq. (jSJ) with 16m c (d)y/9 and by replacing my with 2m c . (In the calculation of Ref. |2j, 
the relative momentum of the c and c that form the charmonium and their binding energy 
were neglected, and so the invariant mass of the cc pair was taken to be 2m r .) In Table II of 

Q 

Ref. H, the cross section for J/ip + J ftp was given as 6.65 ± 3.02 fb, where the uncertainty 
comes only from varying the charm quark mass over the range 1.2 GeV < m c < 1.6 GeV. 
This result includes both the fragmentation diagrams in Fig. [2 and the nonfragmentation 
diagrams in Fig. |21 It does not include QCD radiative corrections or relativistic corrections, 
which are known only for the fragmentation term. As is described in Ref. Q], if the known 
corrections to the fragmentation term in the cross section are applied to the entire cross 
section, then the central value for the J/ijj + J/ip cross section decreases to about 2 fb. This 
value is reasonably close to the result in Table ITU but there are large uncertainties that arise 
from the uncertainty in m c . Furthermore, this procedure suffers from the deficiency that 
the radiative and relativistic corrections to the fragmentation term in the cross section are 
not necessarily valid for the entire cross section. 

A more correct procedure would be to apply the known corrections to the fragmentation 
amplitude to that amplitude alone in the calculation of Ref. |2|. However, this approach 
would still be hampered by a large uncertainty in the fragmentation amplitude that arises 
from the uncertainty in m c . Furthermore, there are higher-order radiative and relativistic 
corrections to the fragmentation amplitude that would not be included in this approach. 
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Both of these drawbacks can be eliminated by writing the fragmentation amplitude in 
terms of gy 1 and my instead of {0%)y and m c . The resulting expression for the fragmenta- 
tion term in the differential cross section is Eq. (j5J), which has no explicit dependence on the 
charm-quark mass. Note that this procedure cannot be used to account correctly for the ra- 
diative and relativistic corrections to the nonfragmentation amplitude. These corrections to 
the nonfragmentation amplitude do not necessarily correspond to those that are contained 
in gv-y. Furthermore, it would be inappropriate to set m c equal to my/2 in the nonfrag- 
mentation amplitude, since the vector-meson mass has little to do with the propagation of a 
charm quark at short distances of order \jE\, eaxa . In the case of production of J/ij; + ij;(2S), 
no consistent choice would be possible because mj^ ^ m^s)- 

We now explain precisely how we calculate improved cross sections for two charmonium 
vector mesons that take into account the nonfragmentation contributions. For each pair of 
vector meson helicities Ai and A2, the cross section can be expressed as the product of a flux 
factor l/(2s), the square of a T-matrix element, and a phase-space factor A 1 ^ 2 (l, r\, ?~2)/(87r). 
We use the physical vector meson masses in the phase-space factor. The T-matrix element 
is the sum of a fragmentation amplitude and a nonfragmentation amplitude. The cross 
section is the sum of a fragmentation contribution, a nonfragmentation contribution, and an 
interference contribution. Our expression for the fragmentation amplitude is proportional to 
9Vn9v2-yi with a coefficient that is a function of the meson masses my 1 and my 2 . Hence, the 
fragmentation contribution to the cross section is given by Eq. Our expression for the 
nonfragmentation amplitude is proportional to {Amy 1 my 2 (0)y 1 (0)y 2 Y^ 2 , with a coefficient 
that is a function of the quark mass m c . The factors (2my i ) 1//2 arise from the relativistic 
normalizations of the meson states. Aside from these normalization factors, the square of 
the nonfragmentation contribution to the T-matrix is identical to that in Ref. |2j. The 
interference contribution to the square of the T-matrix element depends on both the meson 
masses and the quark mass. We calculate the cross sections for each helicity combination and 
then add them. This approach is convenient for taking into account the difference between 
the polarization vector for a meson with mass my i and the polarization vector for a quark 
pair with invariant mass 2m c . The interference and nonfragmentation contributions to the 
cross section are given in Eqs. (|A5|) and (|A6|) in the Appendix. 

Next, let us describe how we estimate the residual theoretical uncertainties. The un- 
certainties are obtained by adding five theoretical uncertainties in quadrature. The only 
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significant uncertainties in the fragmentation amplitudes are those that arise from the elec- 
tronic widths of the vector mesons, which enter through the overall factor of gv-adv-w These 
uncertainties affect the fragmentation and interference terms in the cross sections. The most 
significant uncertainties in the nonfragmentation amplitudes arise from the NRQCD matrix 
elements, QCD radiative corrections, relativistic corrections, and the charm-quark mass. 
These uncertainties affect the interference and nonfragmentation contributions to the cross 
sections. We estimate the error associated with the charm-quark mass by varying m c in the 
interference and nonfragmentation terms in the cross section over the range m c = 1.4 ± 0.2 
GeV. Our estimates of the uncertainties in the NRQCD matrix elements (0\)v are de- 
scribed below. We assume that the QCD radiative corrections to the nonfragmentation 
amplitude are of relative size a s (2m c ) = 0.25. We assume that the relativistic corrections 
are of relative size {(v 2 )v 1 + ( v2 )v 2 )/2- Radiative corrections to the fragmentation ampli- 
tude that involve the exchange of hard gluons between the two mesons are suppressed as 
a^(i?beam) ( m c/-^beam) ~ 0.4% and can be neglected. Corrections that involve the exchange 
of soft gluons between the two mesons are suppressed as (m c t>) 4 /i?^ eam « 0.04% and can 
also be neglected. 

We compute the NRQCD matrix elements (Ci)jm and (Oij^S) from the electronic 
widths of the charmonia by making use of the formula 



which follows from Eq. (Jz3J). Note that Eq. (|25|) differs from the NRQCD factorization 
formula that is given in Ref. Q|: m c has been replaced with my/2 and, consequently, the 
relativistic correction has changed from — |(t> 2 ) v to — \{v 2 )v- The formula (|25[) is less subject 
to uncertainties in the value of (v 2 )v than the standard NRQCD formula and, through its 
dependence on my, resums some corrections of higher order in v. In our calculation, we take 
(v 2 ) J/lP = 0.25 ± 0.09 (Ref. Q), (v 2 )^(2S) = 0.45 ± 0.19 (Ref. Q), and a s {2m c ) « 0.25. 
We obtain 

(Oi)j/^ = 0.482 ± 0.049 GeV 3 , 
(Ci)^ ( 25) = 0.335 ± 0.080 GeV 3 . (26) 

The error bars are obtained by combining in quadrature the uncertainties from the elec- 
tronic widths, the uncertainties from the values of (f 2 )y i; and the estimated errors from 
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TABLE III: Cross sections in units of fb for e + e~ — ► V\ + Vi at -Ebeam = 5.29 GeV for charmonium 
vector mesons, calculated by including both the fragmentation and nonfragmentation amplitudes, 
as is described in the text. The four rows give the fragmentation, interference, nonfragmentation, 
and total contributions to the cross sections. The uncertainties are obtained by combining five 
uncertainties in quadrature, as is described in the text. 



cross section 


J/ij> + J/if> 


J/ij> + </>(2S) 


^(25) +^(25) 


fragmentation 


2.52 ±0.13 


1.81 ±0.06 


0.32 ±0.02 


interference 


-0.98 ±0.48 


-1.09 ±0.60 


-0.30 ±0.19 


nonfragmentation 


0.15 ±0.16 


0.23 ±0.29 


0.09 ±0.14 


total 


1.69 ±0.35 


0.95 ±0.36 


0.11 ±0.09 



uncalculated radiative and relativistic corrections, which we assume to be of relative size a 2 
and (v 2 )y in the rate, respectively. Note that the values of the NRQCD matrix elements in 
Eq. (|26|) are considerably larger than those that were used in Refs. ^. 

Our results for the cross sections for the production of J/ip + J/ip, J/ip + ip(2S), and 
ip(2S) +ip(2S) are given in Table ITTT1 Note that the error bars for the total cross section are 
less than the error bars for the fragmentation, interference, and nonfragmentation contri- 
butions added in quadrature because the interference uncertainties are 100% anticorrelated 
with the fragmentation and nonfragmentation uncertainties. The cross section for J/ip + J/ip 
is 1.69 ± 0.35 fb. The central value is 33% smaller than the value in the fragmentation ap- 
proximation f Table ILTl or the first row of Table ITLTj) . The 33% correction to the fragmentation 
approximation is much larger than the estimate rnl/E^ CSuTa « 7% that was given at the end 
of Sec. IIHI The reason for the large relative size of this correction is that radiative and 
relativistic corrections to the fragmentation contribution to the cross section, which are con- 
tained implicitly in the quantities gy a and my i , decrease the fragmentation contribution by 
about a factor of three from its value without radiative and relativistic corrections. We are 
able to obtain a reasonably accurate result for the J/ip + J/ip total cross section, in spite of 
the large relative errors in the nonfragmentation and interference terms, because the total 
cross section is dominated by the fragmentation term. 
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V. COMPARISON WITH EXPERIMENT AND PREVIOUS THEORETICAL RE- 
SULTS 



The BaBar Collaboration has recently measured the cross sections for production of 



p° + p° and p° + in e + e collisions at energy E heam = 5.29 GeV |l7| : 

a[e + e- -> p° + p°] = 20.7 ± 0.7 stat ± 2.7 syst fb, (27a) 
afe+e- -> p° + 0] = 5.7 ± 0.5 stat ± 0.8 syst fb. (27b) 

The measured cross sections are subject to the cuts | cos^l < 0.8, 0.5 GeV < m p < 1.1 GeV, 
and 1.008 GeV < m < 1.035 GeV. 

In Table HVI we give our predictions for the cross sections that involve light mesons, inte- 
grated over the region | cos#| < 0.8. In Table IVl we give our predictions for the charmonium- 
charmonium cross sections, integrated over the same angular region. The cross sections in 
Tables HVl and IVl are much smaller than those in Table HP and Table ITTT1 respectively, because 
the cut on 9 excludes most of the peak near \x\ = 1. In the last row of Table HV| the m p cut 
is implemented. We can calculate the effect of the cut as follows. The cross section to 
produce a meson and one of the other vector mesons is a constant, to within 1.5%, over 
the range 1.008 GeV < < 1.035 GeV. Hence, we can calculate the fractional correction 
from the cut, with errors of less than 1%, by calculating the integral with respect to m| 
of the line shape over the range 1.008 GeV < < 1.035 GeV and comparing it with the 
integral of the line shape from threshold to infinity. Using a simple Breit-Wigner line shape 
and taking the values of and in Tabled we find that the fraction 0.898 of the integral 
of the line shape is contained in the range 1.008 GeV < < 1.035 GeV. Therefore, in 
order to compare the results in Table IIVI with the BaBar results, one should multiply each 
entry in Table IIVI by a factor 0.898 for each meson in the final state. Specifically, we 
obtain <j[e + e~ — > p° + 0] = 5.04 ± 0.18 fb with the mass cuts on both m p and m^. Our 
predictions for p° + p° and p° + in Table IIVI agree with the BaBar results in Eq. (|2"Tj) to 
within the experimental and theoretical uncertainties. 

Our results in Tables ITT1 and IIVI differ somewhat from those in Ref. 111! ■ Some of the 



differences arise because we are using the 2006 compilation of the Particle Data Group 
0| , rather than the 2004 compilation Q], for the meson masses and widths. However, 
once these differences in the input data are taken into account, discrepancies still remain. 
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TABLE IV: Cross sections in units of fb for e + e~ — ► V\ + Vi at -Ebeam = 5.29 GeV, calculated 
using the fragmentation approximation. The angular cut | cos#| < 0.8, which is used in Ref. 
has been applied. The uncertainties shown are only those that arise from the uncertainties in the 
electronic widths of the vector mesons. The first three rows are calculated in the narrow-width 
approximation. The last two rows are calculated by taking into account the nonzero width of the 
p meson, as is described in the text. In the last row, the cut 0.5 GeV < m p < 1.1 GeV, which is 
used in Ref. 17], has been applied. 



Vx \ V 2 


P° 







J/il> 


V>(2S) 


P° 


23.30 ± 0.80 


3.93 ±0.15 


6.43 ± 0.23 


10.92 ±0.33 


4.44 ±0.13 






0.17 ±0.01 


0.54 ±0.02 


0.92 ± 0.04 


0.37 ±0.02 








0.44 ±0.03 


1.50 ±0.06 


0.61 ±0.02 


p° (no mass cut) 


20.50 ±0.71 


3.68 ±0.14 


6.03 ± 0.22 


10.24 ±0.31 


4.17 ±0.12 


p° (mass cut) 


17.71 ±0.61 


3.42 ±0.13 


5.61 ±0.20 


9.52 ±0.29 


3.87 ±0.11 



TABLE V: Cross sections in units of fb for e + e~ — ► V± + V% at i^beam = 5.29 GeV for charmonium 
vector mesons, calculated by including both the fragmentation and nonfragmentation amplitudes, 
as is described in the text. The angular cut | cos#| < 0.8 has been applied. The four rows give 
the fragmentation, interference, nonfragmentation, and total contributions to the cross sections. 
The uncertainties are obtained by combining five uncertainties in quadrature, as is described in 
the text. 



cross section 


J/ij> + J/V 


J/il> + ifi(2S) 


ii>(2S) ±V(2S) 


fragmentation 


1.20 ±0.06 


0.94 ± 0.03 


0.18 ±0.01 


interference 


-0.72 ±0.36 


-0.81 ±0.45 


-0.22 ±0.15 


nonfragmentation 


0.13 ±0.14 


0.20 ±0.25 


0.08 ±0.12 


total 


0.60 ±0.24 


0.33 ±0.24 


0.04 ±0.06 



The largest discrepancy is between the results for p° + ip{2S) production with nonzero p- 
meson width and |cos^| < 0.8. This discrepancy is about 27%. Other discrepancies for 
cross sections computed with the cut | cos#| < 0.8 are approximately 4% or less. For cross 
sections computed with | cos#| < 1.0, the discrepancies large as 10%. 
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The production cross sections for two vector-meson charmonium states were also calcu- 
lated in the fragmentation and narrow- width approximations in Ref. [[IS]. In this paper, 
somewhat different values for the electronic widths of the charmonium states were used, and 
the resulting production cross sections are lower than ours by about 10-30%. 

The Belle Collaboration has set an upper limit on the cross section for J/ip + J/ip in 
e + e~ collisions at energy -Ebeam ~ 5.29 GeV [20] ]: 

a[e + e- -> J/ip + J/ip] x B >2 [J/ip] < 9.1 fb (90% C.L.), (28a) 
a[e + e- -> J/ifj + ip(2S)] x B >2 [ip(2S)] < 5.2 fb (90% C.L.), (28b) 

where £>>2[V] is the branching fraction of V into final states with more than two charged 
tracks. These upper limits are compatible with the predictions in Table IHII By adding up 
exclusive branching fractions for J/ip decays [9], one can show that the branching fraction 
in Eq. ()28a|) satisfies 13% < B >2 [J/ip] < 80%. It should be possible to measure B >2 [J/ip] 
and B >2 [ip(2S)) at CLEOc or BESIII. The upper limit in Eq. ()28a|) was obtained with a 
data sample of 155 fb -1 at or near the T(4S'). The combined data samples of the Belle 
and BaBar experiments now exceed 1000 fb -1 . Our prediction for the cross section for 
e + e~ — > J /ip + J/ip indicates that there is a possibility that this process can be observed at 
the B factories. 



APPENDIX A: INTERFERENCE AND NONFRAGMENTATION CONTRIBU- 
TIONS TO THE CROSS SECTIONS 

In this appendix, we give the interference and nonfragmentation contributions to the 
differential cross section for e + e~ — > Vi(Ai) + V 2 (X 2 ), where V\ and V 2 are charmonium 
vector mesons, and Ai and A2 are their helicities. The differential cross section for the 
fragmentation contribution can be found in Eq. (jSJ). The interference contribution to the 
differential cross section is 

dahlt \V(\ ) l V(\ )] 10247r 3 a y 7 ^ 27 (4e> yi m V2 (O) yi (O) y2 ) 1 / 2 A 1 / 2 (l,^,r 2 2 ) 
dx [ 1{ l) 2[ 2)1 " 3s 5 r 2 r 2 r 2 (l - r\ - r 2 ) 2 [l - (1 - A)x 2 ] 

x^ A2 (n,r 2 ,r,x), (Al) 

where r = 4m c /y / i. (Note that r± and r 2 reduce to r/2 in the nonrelativistic limit.) In 
Eq. (jAl|) . the numerator factor in parentheses corresponds to the expression for gy^gy^ at 
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leading order in a s and v [Eq. Similar factors appear in subsequent equations for cross 

sections in this Appendix. The functions (ri, r 2 , r, x) are given by 

Fi n * ±1 (ri,r 2 ,r,x) = rV(l - x 2 )[(l - r{ - r\) - A(l, r 2 , r 2 2 )], (A2a) 

fS T i(n,r 2 ,r,x) = (l-r 2 -r 2 )(l-x 4 ), (A2b) 

^±i,o(n,r- 2 ,r,2;) = rr 2 [(l - r\ - r 2 )(l - x 2 )(l - 2x 2 ) + 4r 2 x 4 ], (A2c) 

^Si(n,r 2 ,r,x) = rri[(l-r 2 -r 2 )(l-a; 2 )(l-2x 2 )+4r 2 x 4 ], (A2d) 

^o3(n, r 2 , r, x) = 4nr 2 (l + r 2 )x 2 (l - x 2 ). (A2e) 

The nonfragmentation contribution to the differential cross section is 

d* ni U/n ^ T/ „ v, _ 81927r 3 a 4 (4e 4 m yi m y2 (0) yi (Q) y JA 1 / 2 (l,r 2 ,r 2 )F° f iA2 (r,x) 

— ^(AO + VatAa)]- ^ , (A3) 

where 

F± lt±1 (r, x) = 2r 4 x 2 (l - a; 2 ), (A4a) 

F£{ Tl (r,x) = 2(l-x 4 ), (A4b) 

Ff 10 (r, x) = F f ±1 (r, x) = r 2 (l - 3x 2 + 4x 4 ), (A4c) 

F M (r, x) = 2(1 + r 2 ) V(l - x 2 ). (A4d) 

After summing over helicity states, one finds that the interference and nonfragmentation 
contributions to the differential cross section become 

da^ 2048<K 3 a 4 g Vll g V2 ,(4e 4 c m Vl m V2 (0) Vl (0) V2 ) 1/2 X 1/2 {1, r 2 , r 2 ) 



dx v n ' " 2 ' c/ 3s 5 r 2 r 2 r 2 (l - r 2 -r 2 ) 2 [l - (1 - A)x 

x 



and 



|4rrir 2 (ri + r 2 ) + (1 - x 2 ) [r 2 (r x + r 2 ) 2 (l - r 1 + r 2 )(l + ri — r 2 ) 
— r(ri + r 2 )(l - r 2 + 8rir 2 - rf) + 2(1 - r 2 + r x r 2 - r 2 )] 
-(1 - x 2 ) 2 [l - r(n + r 2 )] 2 (l - n + r 2 )(l + n - r 2 ) I (A5) 



da 11 ^ , 163 847r 3 a 4 (4e 4 my 1 m 1 / 2 (0)y 1 (0)y 2 )A 1 / 2 (l,r 2 r. 

-(m Vl ,mv 2 ,m, 



dx 11 2 ' 9s 5 r 4 

x [4r 2 + (1 - x 2 )(3r 4 - 8r 2 + 5) - 3(1 - x 2 ) 2 (l - r 2 ) 2 ] . (A6) 
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